Re-entrant ordering of solute in a colloidal suspension during solvent evaporation 
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We study the phenomenon of self-assembly of silica micro-spheres on a glass plate during evapo- 
ration of the solvent from a colloidal suspension. Our experiments unveil an interesting competition 
between ordering and compaction in a strongly driven, out of equilibrium system arising from a 
slowing down of dynamics due to an impending glass transition. A suitable choice of experimental 
conditions minimizing the influence of many other competing phenomena that usually complicate 
probing of this underlying physics is crucial for our study. A re-entrant behavior in the order-disorder 
phase diagram as a function of particle density and drying time is established and the results are 
explained with the help of simulations and phenomenological theory. 

PACS numbers: Valid PACS appear here 



Self-assembly of colloidal particles into well-ordered ar- 
rays is interesting both due to the diverse technologi- 
cal possibilities [TJ |2] and because it involves ordering in 
systems driven far from equilibrium [3 . Ordered, two- 
dimensional particle arrays are often produced by evapo- 
rating a drop of colloidal suspension on a substrate [4j [5]. 
The physics of ordering upon evaporation is poorly un- 
derstood because of the complex interaction of many 
forces involved viz. lateral capillary forces, convection 
and floatation, as well as colloid-colloid interactions, and 
contact line pinning [6-11 !. In this Letter, we report an 
experimental realization that enables us to study the un- 
derlying interesting physics, including reentrant behavior 
of such driven systems as a function of just two easily con- 
trollable parameters, namely the initial particle density, 
po and the drying time, of the solvent, the ordering 
mechanism being essentially reduced to a competition 
between drying and slowing down of dynamics due to an 
impending glass transition[12 . 

It is commonly believed that drying assisted assembly 
of colloidal particles proceeds by a two-stage mechanism. 
In the pinned state, the contact line at the periphery of 
the drop is pinned by the substrate, the drop area is fixed 
and a flow, which advects particles to the periphery, is 
set up to replenish the liquid lost due to evaporation [7J[8] 
which occurs primarily at the periphery U due to the in- 
creased local curvature. The build-up of particles at the 
periphery causes self-assembly due to attractive capil- 
lary forces, that has been investigated in the past [8| [TO] . 
However, the contact line is eventually depinned and dry- 
ing recommences when the drop area decreases [TT] . The 
phenomenon of self assembly in this regime has not been 
investigated so far. The lifetime of the droplet in the two 
regimes and, consequently, their relative importance de- 
pend crucially on the contact angle 6 C between the drop 
and the substrate. For small C , the lifetime of the pinned 
stage is small. We have used this fact to our advantage by 
decreasing C below 25° by preparing highly hydrophilic 



substrates, enabling the suspension to dry predominantly 
under de-pinned conditions. Additionally, because of the 
low value of C , the evaporation rate from the periph- 
ery is weak making the outward flow rate of the liquid 
and hence the centrifugal particle flux negligible. In or- 
der to ascertain the irrelevance of hydrodynamic effects, 
such as Marongoni flows [6 , in determining the final con- 
figuration of the silica particles, we estimate the Peclet 
number, P n = r p v p /D c1 to be less than 0.02 in all cases, 
where, r p is the radius of the particle moving with a ve- 
locity of v p and D c is the diffusion constant for all and 
po that we have experimentally investigated. 

Colloidal suspension of silica micro spheres of radius 
0.554(±0.023) /im were obtained from Bang's Labora- 
tory, USA. In order to prepare substrate with a highly 
hydrophilic nature glass slides (cover slips) were heated 
at 70° C for 6 hours in a piranha solution (3 : 1 mixture of 
cone. H2SO4 and H2O2). Then they were cooled, washed 
several times with de-ionized water and dried. In each 
set of experiments freshly prepared glass slides were used. 
Slides obtained by this method are highly hydrophilic 
with a contact angle less than 25° in all cases. In a typ- 
ical experiment, 60 jA of silica suspension in deionized 
water with different silica concentrations (vol%), po? each 
time was drop-cast on a substrate and dried upto a total 
drying time, td, till all trace of the solvent evaporated. 
Humidity and temperature are the available variables to 
obtain different values of td over a wide range. The tem- 
perature was 25 K for all measurements, except when 
we needed to achieve the shortest td (30 mins.) used in 
our experiments, when the temperature was set at 35 K. 
It is known that treatment of glass slide with piranha 
solution makes the surface of the glass substrate nega- 
tively charged. Zeta potential measurements show that 
the surfaces of these colloidal silica particles are also neg- 
atively charged (Zeta potential -41.3 mV). Therefore, 
the interaction between the substrate and the particles 
is repulsive, reducing random pinning of particles on the 
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FIG. 1: (a)and(b): Real space SEM image of silica sphere, Initial volume fraction = 0.15%, (scale bar = 20 /xm) with a drying 
time of 30 min and a drying time of 14 hours, respectively. Inset shows the higher magnification image (scale bar 2 pm) (c) 
and (d): Corresponding structure factors, (e) and (f): Change in the average domain size D av and bond order parameter 
Tne(pojtd)j respectively, as a function of different silica concentration, p at different drying time td- (g) and (h): Change in 
the average domain size D av and mo(po,td) as a function of td at different po respectively, a typical error bar in term of the 
standard deviation calculated from different runs of the measurements is shown as a bar in each of fig (e) to (h) 



substrate allowing us to study a cleaner limit of the self- 
assembly under a driven condition. The interaction be- 
tween the particles depends not only on their surface 
charges, but also on the ionic strength of the solvent, 
which screens the surface charge thus playing an impor- 
tant role in controlling the interaction between particles 
and hence the ordering process [13]. The screening length 
of the solvent was kept fixed for all the runs, thereby re- 
ducing one more controlling parameter in the problem, 
by fixing the ionic strength of the solvent by the addi- 
tion of a known concentration of KC1 (10 x -5 mols/lt) 
in each case. The random percentage fluctuation of the 
screening length could thus be reduced far bellow what 
can be achieved with de-ionized mili-Q water (resistance: 
18 Mohm). 

Scanning Electron Microscope (SEM) pictures were 
taken of the final structure of the colloidal particles after 
the completion of the drying process using a FEI- Quanta 
200. Two representative images are shown in (Fig.l (a) 
and (b)). To quantify the extent of ordering, coordi- 
nates of every particle in a snap shot of the SEM image 
over a large area (50 pm x 50 pm) were obtained using 
an image digitizer software (Image J). We characterize 
each such configuration by calculating the structure fac- 
tor S(q) = N~ x J2 m n exp(i q • r mn ), with the sum taken 
over all pairs of particles index by m and n and separated 
by distance r mn = r m — r n . S(q) is isotropic for disor- 
dered configurations and shows crystalline peaks for or- 
dered ones (Fig.l (c),(d)). The average domain size, D av 



may be obtained from S(q) using the well known scher- 
rer formula [T4]. The orient ational (Nelson-Halperin) [T5] 
order parameter ra 6 = N~ x z\ Y,k=i Y? n =i ex P(^ 6# fen )|, 
where N is the total number of particles and z the coor- 
dination number, was also computed to characterize the 
final configurations. (Fig.l (e),(f)) show the variation of 
D av and m§ with po for various fixed drying times, td- 
It is clear from the results that for slow drying, there 
is a jump in both D av and m§. Whereas, in case of 
faster drying no such transition is observed in either of 
the two parameters characterizing the extent of ordering, 
with increase in the concentration of the initial suspen- 
sion. Similar results are also obtained when the data are 
plotted against td at different po (Fig.l (g),(h)). 

Guided by the jump, we use threshold value of 5 /am 
and 0.69 for D av and m$ respectively, (marked by hor- 
izontal lines in each of (Fig.l (e)to(h))) to distinguish 
between ordered and disordered phases, allowing us to 
obtain a dynamical phase diagram in the td — po space 
for this whole process (Fig. 2). We note that results pre- 
sented here are not critically dependent on the exact 
choice of the thresholds and a range of values for the 
thresholds leads to same results. (Fig. 2) clearly shows 
that the system is invariably disordered for po < 0.1% 
of silica particle. For po > 0.1%, we find that the sys- 
tem is disordered for small td and is ordered for larger td] 
the critical td to bring about the ordered structure in- 
creases somewhat for a higher po. Thus the data clearly 
establish a disorder to order transition as a function of 
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po across 0.1% for a wide range of drying time td > 4 
hrs. Likewise, there is a disordered to ordered transition 
as a function of the increasing drying time td for po > 
0.1%. Interestingly, present data also show an ordered to 
disordered transition on increasing po for lower values of 
td < 10 hrs. The data points shown in (Fig. 2 ) therefore 
establish a dynamical phase diagram with a re-entrant 
disorder — ^ order — ^ disorder transition in this system, 
the phase boundary going through a minimum showing 
that it is possible to optimize both td and po in order to 
obtain large, ordered, self-assembled arrays. 
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FIG. 2: The experimentally obtained dynamical phase dia- 
gram in the td — po plane. The red asterisk and the green 
open circle represent ordered and disordered configurations, 
respectively. Overlapping asterisks and circles denote points 
at the phase boundary, where independent experimental runs 
yielding contrasting extent of ordering across the threshold, 
indicative of the existence of the phase boundary. Also shown 
are results from our BD simulations: blue circles. We use 4000 
particles for our simulations and a maximum of 625 cells each 
of size 5 to represent the solvent field 0. A time step of 10 -3 is 
used to update both particle positions and <j). Comparing our 
numerical solution with the analytic result for Ising droplets, 
we estimate our accuracy for updates to be 1 in 10 -6 . In 
order to compare our results with experiments, we have mul- 
tiplied the axes poK 2 by 0.2 and jtd by 1.5 to obtain po in 
units of Vol. % of silica and td in hours. The scale factor for 
po was obtained from an independent measurement of k for 
the silica particles in solution while td was simply scaled to 
obtain best fit with experiments. 



To understand the effect of td and po on the final struc- 
ture, we model [16] the drying process as a 2d liquid-gas 
phase transition. The local density of the solvent 0(r, t) 
at position r, satisfies [17 , (j) — 5 T '/5(f) with ^[{^j] = 
J cPr a(j) 2 + (j) 4 + c(V(/>) 2 + fJL 3 (f>, where p s is the chem- 
ical potential of the solvent controlling the rate of dry- 
ing and a,c are parameters. The homogeneous part of 
the free energy density has minima at 0o = ±V — 2a/2 
when a < and p s = denoting the liquid (j) = <po 
and vapour <\> = — <po phases of the solvent. The parti- 
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FIG. 3: (a) Data collapse plot of the bond order parameter 
ms/m^ vs td/r for various values of td and po; mg° is the 
saturation value of mo and r is the (po dependent) relax- 
ation time, (b) Plot of the relaxation time r vs poK 2 showing 
a minimum. Open circles represent data from simulations 
as obtained from the scaling plot in (b), where the line is a 
fit with the analytic form (see text).(c) Plot of the average 
coordination number z (open squares, scale to the left) and 
fractional increase in the density Sp/po (filled squares, scale 
to the right) as a function of poK 2 for fixed td — 7 hours. 



cle dynamics is modeled using an over-damped Brow- 
nian dynamics (BD) scheme, r(0)r^ = i^, where the 
force Fi = ^ ~d^r + J^r + Vi consists of three parts 
(1) from the pair interaction Vij = Vq exp(— K.rij)/rij 
(screened Coulomb) between the particles i and j, (2) 
from surface tension forces produced at the contact line 
and (3) Gaussian white noise rji at the ambient temper- 
ature. In order to model the fact that the silica particles 
are immobile without the solvent, we take the mobil- 
ity T = (tanh(70) + l)/2, which vanishes in the limit 
(j) — >• —0o El]- The prefactor, Vb, sets the energy and 
k the screening length the scale for distances. The con- 
trol parameters are the initial density of the particles in 
the drop, po, and the drying rate which is controlled by 
varying p s . 

For a given po, we prepare an initial circular drop 
(0(r,£ = 0) = </>o for |r| = R and — 0o otherwise) con- 
taining randomly dispersed particles. As the drop evap- 
orates, the shrinking circular contact line (</>(r, t) = 0) 
draws the particles together increasing the particle con- 
centration. When (j) = — everywhere all dynamics 
ceases (T(— <po) = 0). In order to make a direct compar- 
ison with experimental results we calculate ttiq, taking 
ttiq = 0.69 as the threshold for ordering. The resulting 
phase line is shown in (Fig. 2). Considering the simple 
model that we have used, the agreement between the 
calculated phase line and the one suggested by the ex- 
perimental results, including re-entrance, is remarkable. 

The average order parameter scales as me(p = 
7^6° /(W r ) 5 where rag , the saturation value and the 
relaxation time r are functions of po alone, as illustrate 
by the collapse of the simulated rae/rag as a function 
of (td/r) for several po values shown in (Fig. 3(a)). It 
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is clear that r itself goes through a minimum (see Fig. 
3(b)). The reentrant behavior as seen in the dynami- 
cal phase diagram is thus coded into the shape of r(po). 
Indeed, the dynamical phase diagram is a plot of the con- 
tour rriQ = constant (= 0.69 in our case) being given by 
td = r(po) x / _1 (m 6 ). 

There are two timescales associated with the relax- 
ation of the order parameter r = t\ + r^. The time, 
Ti, required for the concentration to increase from the 
initial value po to the level pf the freezing density, can 
be easily estimated [18], being simply given by the time 
taken by an Ising droplet to shrink by a fixed amount viz. 

-i 3/2 

r\ = Aip + A2p + As where Ai depend on the pa- 
rameters c and n s as well as pf. Note that T\ diverges as 
po — >• 0. The time required for order to develop, T2, is the 
time needed by the system for exploring the space of all 
configurations. A Vogel-Fulcher like empirical form [12] 
for the relaxation time, viz. T2 = Toexp[C/(pj — po)] 
with To, C and pj as free parameters, fits our simulation 
data rather well (Fig. 3(b)). For small po, fluctuations al- 
low the system ample time to explore configuration space 
during the drying process, though one requires a large 
drying time (td(^i) m order to raise the density of the 
solution to the value needed to initiate ordering. On the 
other hand if po is high, rapid compaction causes the par- 
ticles to jam, setting up a network of force chains [T9H2T] 
which resists any further increase of density. Ordering in 
this limit is dominated by T2. Slowing down of the dy- 
namics is caused by divergences of t\ and T2, alternately 
at the low and high density limit, causing the reentrant 
transition observed in these experiments. Deep inside 
the high density region an additional complication, viz. 
the formation of multiple layers commences in the exper- 
imental system cutting off the divergence of T2. 
In this Letter, we have shown that drying induced order- 
ing in carefully prepared samples can be used to study 
in detail the competition between the thermal and ather- 
mal processes of ordering and compaction. Minimizing 
additional complications that often arise in most experi- 
mental systems, eg. multilayer formation, hydrodynamic 
effects, pinning and poly-dispersity, which can be investi- 
gated separately in the future, our experiments, nonethe- 
less, show a rich phase diagram, including a clear sign of 
reentrant behavior as a consequence of this competition 
at the simplest level. We explain these phenomena in 
this system driven far from equilibrium as arising from 
a slowing down of dynamics due to an impending glass 
transition in terms of simple physical arguments based 
on detailed Brownian dynamics simulations. 
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